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Introduction 


Let  an  analytic  function  f  of  a  single  complex  variable  be  defined 
in  a  neighborhood  of  a  point  Zq  by  means  of  its  Taylor  series  at  z^. 

In  this  paper,  we  wish  to  discuss  a  constructive  method  for  the  solution 
of  the  following  problem:  Suppose  it  is  known  that  the  function  f  can 
be  continued  analytically  into  a  domain  R  of  the  complex  plane.  It  is 
desired  to  compute  the  value  of  f  at  an  arbitrary  point  b  e  R, 

By  a  constructive  method  (or  "algorithm")  for  solving  a  problem 
in  numerical  analysis  we  mean  a  sequence  of  rational  functions  of  the 
data  of  the  problem  which  converge  to  the  desired  mathematical  object. 

In  a  number  of  special  situations  classical  analysis  offers  a 
variety  of  formulas  and  methods  for  solving  the  problem  posed  above, 
such  as  the  Schwarz  reflection  principle,  or  the  methods  of  Mittag- 
Leffler  and  Borel  (see  Bieberbach  [1955]).  The  problem  is  easily  solved 
if  the  function  f  is  known  to  satisfy  a  differential  or  other  functional 
equation.  In  some  cases  the  problem  may  also  be  solved  by  classical 
methods  of  summability  theory  (see  Hardy  [1948]).  In  a  more  numerical 
vein,  Euler's  series  transformation  sometimes  proves  to  be  an  effective 
means  of  continuation  (van  Wijngaarden  [1953]).  Recently,  linear  pro¬ 
gramming  techniques  have  been  proposed  for  the  solution  of  a  special 
continuation  problem  (Douglas  [i960];  Douglas  and  Gallie  [1959]). 

The  above  methods  apply  only  in  special  (sometimes  extremely  so) 
situations.  However,  the  following  technique,  due  to  Weifers trass,  is 
generally  applicable  (for  a  recent  exposition,  see  Behnke  and  Sommer 
[1955] ,  pp.  166-171)  .  We  join  the  points  z^  and  b  by  a  rectifiable 
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arc  Y  lying  in  R  .  If  6  >  0  denotes  the  distemce  of  y  from 
the  boundary  of  R  ,  we  select  N  points  z^,z^, , . . fZ^  =  b  on  y 
such  that  |zj^  -  I  <  k  =  1,2,...,N.  The  Taylor  series  of  f 
at  Zq,  having  a  radius  of  convergence  >6,  enables  us  to  calculate 
the  values  of  f  and  of  all  derivatives  of  f  at  z^.  Thus  the  Taylor 
series  at  z^  is  also  known,  and  can  be  used  to  calculate  the  coefficients 
of  the  Taylor  series  at  the  point  z^.  Proceeding  in  this  manner,  we 
finally  obtain  the  Taylor  series  at  the  point  z^^,  and  thus  a  fortiori 
the  solution  of  our  continuation  problem. 

Clearly,  the  above  method  is  not  constructive  in  the  sense  indicated 
earlier.  The  data  in  this  case  are  the  coefficients  of  the  Taylor  series 
at  z  =  a.  While  it  is  tine  that  the  partial  sums  of  the  first  Taylor 
series,  and  its  derivatives,  are  rational  functions  of  these  coefficients, 
their  limits  generally  are  not.  Yet  these  limits  are  required  to  continue 
the  process .  In  what  follows  we  shall  describe  a  modification  of  the 
above  Weierstrass  ian  method  of  continuation  which  transforms  it  into  a 
constructive  process. 

2.  Matrix  formulation  of  Weierstrassian  analytic  continuation. 
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Let  a^  +  a^(z  -  Zq)  +  a^Cz  -  z^)  +•••  be  the  power  series 

defining  f  near  Zq,  and  let 

f(z)  =  tiQ  +  b^(z  -  z^)  +  b^Cz  -  z^)^  +••• 

near  z^.  If  i z^  -  Zq |  is  smaller  than  the  radius  of  convergence 
of  the  series  at  z^,  we  have 

•>0  =  *  '■1*2  + 
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3 


where  “  ^0’  more  generally  for  k  =  0,1,2,.. 


,  ,  V  ..  I  -..y  .  ,  ,2  (k  +  2)  I 

k  ^1  II  ^k+1  ^1  21  ^k+2 


=  (S)“k"('‘r)'‘l\+l^(‘'2')''2 

(:) 


^k+2  ■^* 


where 


denotes  a  binomial  coefficient. 


We  can  render  these  formulas  more  lucid  by  the  use  of  matrix 
notation.  We  denote  by  A(z)  the  infinite  column  vector  whose  components 
are  the  coefficients  of  the  Taylor  series  at  the  point  z.  For  instance, 


We  now  define  an  infinite  upper  triangular  matrix  M(z) 
M(z)  —  m,n  —  0,1,2,..., 


as  follows : 


where 


n  >  m 


n  <  m  . 
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Written  out  in  full. 


The  relation  (l)  now  can  be  combined  into  the  single  relation 

(2)  A(z^)  =  M(6^)  A(zq)  . 

In  a  similar  manner  we  find,  if  \\  ~  I  less  than  the  radius  of 
convergence  of  the  Taylor  series  at  \ 

(3)  A(zj^)  =M(6j^)  A(zj^_^),  k  =  l,2,...,N, 

where  1'  solution  of  our  continuation  problem  is 

given  by 

(4)  A(z^)  =M(6j^)  •••  M(b^)  A(Zq), 

with  the  understanding  that  the  product  is  to  be  formed  proceeding  from 
the  right  to  the  left. 

The  non-associativity  of  the  product  in  (4)  is  shown  by  the  following 
remark;  It  is  easily  seen  that  the  "continuation  matrices"  M(z)  satisfy 
the  addition  theorem 

(5)  M(z')  M(z")  =  M(z'  +  z") 

for  arbitrary  complex  z'  and  z" .  Since  1  '*'***’^  ^1  ~  ~  ^0* 
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the  product  on  the  right  of  (4)  thus  Is  formally  equivalent  to 
M(Zjj  -  Zq)  A(zq).  However,  this  expression  is  without  analytic 
content,  since  it  is  tantamount  to  evaluating  the  Taylor  series  at 
Zq  for  z  =  Zjj,  where  it  diverges  in  general. 

3 •  Definition  of  the  continuation  algorithm 


We  denote  by  finite  matrix  consisting  of  the  first 

m  rows  and  n  colvimns  of  the  matrix  M(z)  .  Similarly,  we  denote  by 
A^(z)  the  column  vector  comprising  the  first  n  elements  of  the  vector 


A(z)  .  By  the  symbol  A^  we  denote  vectors  intended  to  approximate 

A  (  z,  )  . 
n^  k'' 


At  first  sight,  it  would  seem  reasonable  to  transform  the  Weier- 
atraasian  method  of  continuation  into  an  algorithm  by  forming  the 
vectors 


A^^^  =  M  (6  ) 
n  n,n'  N 


M  (V  n) 

n,n'  N-1' 


M  (5-,) 

n,n^  1^ 


V^o’- 


where  n  =  1,2,....  However,  since  products  of  finite  matrices  are 
always  associative,  and  since  the  addition  theorem  (5)  also  holds  for 
the  finite  segments  this  amounts  to  nothing  more  than  to  forming 


=  M 

n,n 
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which  is  precisely  what  we  obtain  by  substituting  z^^  -  z^  into  the 
power  series  defining  f  at  z  =  z^.  It  thus  is  clear  that  if  we  wish 
to  obtain  a  convergent  algorithm,  the  finite  segments  of  the  matrices 
M(6j^)  in  (4)  must  be  chosen  in  a  more  sophisticated  manner. 
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It  is  fairly  obvious  that  if  convergence  is  to  be  assured,  the 
vectors  A(Zj^)  must  be  approximated  particularly  well  if  k  is  small. 
(The  truncation  error  is  propagated  forward  but  not  backward.)  We  can 
achieve  this  by  replacing  the  matrices  t>y  rectangular  segments 

that  have  q  times  as  many  columns  as  rows ,  where  q  is  an  integer 
>  1.  We  thus  are  led  to  defining  the  approximating  vectors 
the  following  manner:  For  n  =  1,2,...,  let 


(6)  =  M  (&,.)  M  ^(6-,  J 

^  '  n  n,nq^  N'  2'  N-1 

nq,nq 


nq  nq 


These  vectors  evidently  depend  on  the  "magnification  ratio"  q.  For 
each  n,  they  can  be  built  up  recursively  as  follows: 


a'°1=a 


(7) 


.N 


nq  nq 

=M  (6)A^^‘^^ 

"  N-k  “  N-k  N-k+l''  V  ^  N-k+1’ 

nq  nq  ,nq  ■  nq 


k  1,2,...,N. 


4.  Convergence  of  the  algorithm. 

Let  rj^  be  the  radius  of  convergence  of  the  Taylor  series  repre¬ 
senting  f  in  a  neighborhood  of  the  point  For  Weierstrassisin 

analytic  continuation  it  is  necessary  that  1 6^^  ]  =  I  I  1’ 

k  =  1,2,...,N.  Let  this  condition  be  met,  and  let  r^^  be  any  number 
satisfying  ^  ^  1’  Cauchy's  theorem  there  exists  a  constant 

C,  such  that 
k 

n  iC 

|f^  ^ (z^  I  <  ^  ,  n  =  0,1,2,.... 

^k 


Setting  Pj^  =  this  may  be  written  in  the  form 

I  ^  °k’  n  =  0,l,2,.... 

k 

Only  a  slightly  stronger  condition  is  required  for  the  convergence  of 
the  algorithm  (7) . 


THEOREM.  Let  there  exist  constants  C-^  and  Pj^,  0<Pj^<l, 
such  that 

(8)  |f(^^(z)  I  <  n!(-|^)  0^,  n  =  0,1,2,..., 

for  all  points  z  on  the  straight  line  segment  .joining  z.^  ^  and 
z^  (k  =  1,2,...,N).  Then  there  exists  a  number  such  that  for 

all  q  >  q^  the  elements  of  the  vectors  converge  to  the  corres¬ 

ponding  elements  of  A(Zj^)  ^  n 

The  question  as  to  the  infimum  of  all  q^  for  which  the  theorem 
is  true  is  left  open.  An  upper  bormd  for  this  infimum,  however,  will 
emerge  from  the  proof. 


5 .  Proof  of  the  Theorem. 


jetting 


=  A  (z  )  -  A^^^ 
n  n  k  n 


(k) 

for  k  =  1,2,.-.,N  and  for  all  values  of  n  for  which  A^  is  defined, 

our  aim  is  to  show  that  R^^^-e  0  as  n  oo.  Qur  first  goal  is  a  recur- 

n 

(k) 

rence  relation  for  R^  .  Such  a  relation  is 

relation  for  the  vectors  A  (z,). 

n  ic 


obtained  via  a  recurrence 
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LEMMA  1.  For  k  =  1,2, ...yN)  n  =  1,2, 

\ 


. ,  let 


dO')  =  p”"! 

n 


/  ("o’) 

(t)  ‘i’ 


Then 


(9)  A  (zj  =  M  (ftjA  (z.  J  + 

'  '  n  k'  n,nq'  k'  nq'  k-1'  k  n  * 

where  denotes  a  diagonal  matrix  whose  elements  are  bounded 

Proof.  By  Taylor's  theorem, 

ST  ^  "  ST  ^  ^  JilnT  ^  ^ 

1  ^(nq-l)^_  s.qn-m-l 

•••+iIT5r^rF^nTT^  (Vl^\ 

\ 

+  -T7 - - - Trr  5 

ml(nq  -  m  -  1)  I  ^  k  ' 


“k-1 


m  =  0,1,2, ... ,n  -  1.  It  follows  that 

A  ( z,  )  =  M  ( 5>i  )  A  ( z,  , ) 
n'  k'  n,nq  k  nq'  k-1 


,nq-l 


+  5  f^^'i^t) 

^k-1 


0!(nq  -  1) 1 

1  t.  j.\nq-2 

II  (nq-  2)1 


\ 


\  Cn  -  1)  l\nq  -  n)  I  / 


(t)dt. 


dt. 


Formula  (9)  now  is  obtained  by  estimating  the  Integral  by  means  of  (8) . 
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Subtracting  (7)  from  (9)  we  now  find 

=  M  (6.)  +  @ 

n  n,nq'  k'  nq  k  n  ’ 

where  n  runs  through  all  multiples  of  q^  A  standard  argument 

In  the  theory  of  difference  equations  now  shows  that 

N 
Z 

k=l  nq‘ 

where  I,  and  for  k  =  1,2,...,N  -  1 


r(N)  ^  I  jy(k)  @  M 

n  ,  ,  '  'n  k  N-k* 


*^n  ^  \,nq^^N^  ^  2^  ^N-1^  ^  N-k-1  N-k^^k+1^* 

'  nq>nq  nq  ,nq 

A 

We  shall  denote  by  A  the  vector  or  matrix  whose  elements  are 
the  absolute  values  of  the  elements  of  A,  and  use  a  notation  such 
as  A  <  B  to  Indicate  the  corresponding  inequality  between  corres¬ 
ponding  elements.  We  then  have,  in  view  of  Lemma  1, 


(10) 


r(n)  <  V  c 


where  ■nl"'  =  I,  and  for  k  =  1,2,...,N  -  1 

=  M  (5j  •••  M  ,,  ,  ^  ,  (6, 

'  'n  n,nq'  N'  N-k-1  N-k'  k+1'^ 

nq  ,nq 

The  resulting  products  of  matrices  in  (lO)  can  be  further  estimated  by 
a  reduction  fonnula  based  on  the  following  result. 


LEMMA  2.  Let 


(11)  <ik  >  1  + 

k  =  1,2,  ...,N  -  1.  There  exist  constants  not  depending  on  q  or 
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n  such  that  for  all  1  and  all  sufficiently  large  n 


where 


(13) 


(q  -  l)'^* 


q-l  ^k^k+l  I  **k+l  '^k  f  • 


-li 


A  value  for  will  b*  given  below. 

Proof.  For  brevity  we  set  |6^[  -  a,  Denoting  by  c 


the  element  of  the  vector 

i,nq<V  °nq’- 


we  have 

(14) 

where 


2 

q  n  -m  „ 

=„  =  Pj  »  3, 


S  . 


nq-lnn 


It  turns  out  that  the  dominant  term  in  the  sum  S  is  the  last  term. 
Factoring  out  this  term,  we  have 


S  -  1  ^  nq  -  m  -  1  a  ,  (nq  ^  m  -  1)  (nq  -  m  -  2) 

^  nq^  -  nq  +  2  ^  (nq^  -  nq  +  2)  (nq^  -  nq  t  3)  ^ 


where 


where 


Since  q  >  2,  m  <  n,  we  have 

The  element  of  being  given  by 


ve  thus  have  from  (14) > (15) , (16)  and  (17) 
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/ 


Estimating  the  factorials  by  means  of  Stirling's  formula,  we  have 


(  nf )  = 


(nq)  l(n(q'^  -  q))  I  •J2r(q 


1 _ r 

1  -  l)n  L(q  - 


and  thus,  in  view  of  the  presence  of  sfn  in  the  denominator, 

,q  Inq 


(19) 


(nq^)  <  r — a! — 


for  all  large  n.  It  thus  follows  from  (l8)  and  (19)  that 


m 


nq 


K,  d  , 
k  m’ 


which  is  the  desired  result. 


It  follows  from  (12)  that 

fr(k)  -(k)  .  nq^“^  fr(k+l)  -(k+l) 

'  'n  °  ^N-k  ^  ^k  I  In  °  N-k-l» 

nq  nq 


and  thus,  using  induction, 

I  In  ^  N-k  ^k+1  ^N-1  \  h+1  ^N-1  n  * 


nq 


:(N) 


We  shall  estimate  by  a  bound  very  similar  to  that  of  Lemma  2. 


LEMMA  3 .  Let 


(21)  qjj>l+|6j^i. 

There  exists  a  constant  independent  of  n  or^  q,  such  that 

A  ( n 

the  elements  of  '  are  bounded  by  for  all  Q  >  and 
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all  large  n,  where 


(22) 


2n  = 


;(N) 


Proof.  We  bound  the  elements  of  the  vector  D'  '  by  the  Siam  of 
-  n 

the  absolute  values  of  all  elements.  This  sum  can  be  formed  by  multi¬ 
plying  by  the  matrix  M,  (l) .  An  application  of  the  method  of 

n  j.  ^ii 

proof  of  Lemma  2  now  yields  the  bound  of  Lemma  3>  where 


%  - 1  -  I'‘nI 


In  view  of  (lO) ,  (20),  and  Lemma  3  we  find  that  every  element  of 


(n) 

R  is  bounded  by 
n 


(23) 


N 

2  c'  ; 
c=l  L 


N-k  ^N-k-1 

1 

c  ^k+l 


I" 

nJ  ’ 


where 

=  ^k  ^k  \il 

and  where  the  numbers  Z-^  are  defined  by  (13)  and  (22)  .  The  proof  of 
the  theorem  now  is  an  easy  consequence  of  the  following  fact: 


LEMMA  4.  For  k  =  1,2, ...,N, 

(24)  lim  \  =  0. 

q-^oo 

Proof.  If  we  set  ^  ^k+1  ^  have 


\  - 1'  Pk  •  =k 
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for  k  =  where 


°k  Pk+1 1  ^k+1 1 


Relation  (24)  now  follows  trivially  from  the  well-known  result 
lim  (1  -  =  e 

q->*oo  ^ 


and  from  the  fact  that,  by  virtue  of  jp^^l  <  1, 


liir  (q  -  1)  pj  =  0. 

q-t-oo 


In  view  of  Lemma  4  there  certainly  exists  q*  such  that  |Zj^|  <  1 
for  q  >  q^^,  k  =  1,2, ..  .,N.  If  q  >  qg,  where  q^  =  max(q»,q^,q2, . . .  ,qjj)  , 
then  all  elements  of  bounded  by  the  expression  ( 23)  ,  which 

tends  to  zero  as  n  This  completes  the  proof  of  the  theorem. 


6.  Choice  of  the  enlargement  ratio 


Once  the  continuation  points  . ,z^,  and  with  them  the 

numbers  Pj^  and  6^^  (k  =  1,2,...,N),  are  fixed,  the  enlargement  ratio 

q  has  to  satisfy  the  following  conditions:  In  order  to  satisfy  (ll) 

and  ( 21)  we  must  have 

6, 

(A)  q  >  1  +  \r^\y  ^  =  1,2, ...,N 

k+1 

(  =  l)  ,  and  in  order  to  make  |Zj^I  <  1,  we  must  require 


(q  -  1) 


q-1  ^k  k 


where  Cj^  is  given  by  (25),  or  approximately 


(26)  (q  -  1)p2  <  k  =  l,2,...,N. 


Condition  (26)  shows  that  imless  is  small,  q  may  have  to  be 

taken  quite  large.  Small  values  of  Pj^,  on  the  other  hand,  force  us 
to  make  N,  the  number  of  continuation  steps,  large.  Thus  in  any  case 
the  continuation  matrices  are  likely  to  be  large.  As  an  example, 

we  consider  the  problem  of  continuing  the  function  log  z  around  the 
unit  circle,  starting  from  z  =  1.  Choosing  the  continuation  points 


z 


k 


iktp  ,  n 

e  ,  where  cp  =  g-  , 


we  find 


Pj^  =  p  =  2  tan  I  =  0.535898. 


Since  all  1 |  are  the  same  in  the  present  case,  condition  (a)  is 
satisfied  for  q  >  2,  and  condition  (b)  simplifies  to 

(q  -  l)p'^”^  <  e“^. 


This  is  satisfied  for  q  >  5,  but  not  for  q  =  4-  Since  N  =  12  in 

the  present  problem,  even  the  very  first  matrix  used  in  the  algorithm  (7) 
12 

has  5  =  244,140,625  columns. 


7.  Numerical  examples 

The  algorithm  (7)  has  been  carried  out  niimerically  in  some  very 
simple  examples.  The  author  is  indebted  to  Mr.  Thomas  A.  Bray  of  the 
Boeing  Scientific  Research  Laboratories  for  his  expert  assistance  in 
the  planning  of  these  computations.  We  wish  to  report  briefly  on  the 
results  of  two  such  computations. 
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(l)  The  power  series 

f(z)  =  1  -  z  +  +•• • 

has  been  continued  using  the  points  Zj^  =  0.35k,  k  =  0,1, 2, 3.  An 

enlargement  ratio  q  =  3  was  used.  The  following  table  lists  the 

(k) 

first  components  of  the  vectors  A'  ' ,  and  of  the  exact  vectors  A(z,). 

XI  iC 


^k 

n 

1 

2 

3 

1 

0.740740741 

0.588239567 

0.492237049 

2 

0.740740741 

0.588235294 

0.487767430 

3 

0.740740741 

0.588235294 

0.487805203 

4 

0.740740741 

0.588235294 

0.487804475 

A(zJ 

0.740740741 

0.588235294 

0.487804878 

The  values  in  the  last  coltimn  approximate  the  function  at  a  point  where 
the  original  power  series  does  not  converge. 

(Il)  The  power  series 

oo  I  ,Nn+l  ^  z  - 

f(z)  =  2  —  z  =  5  log(l  +  t)dt 

n=l  0 


was  continued  using  the  points  =  0.3k,  k  =  0,1, 2, 3, 4.  Again  working 

*  \ 

with  q  =  3,  the  following  values  of  the  first  components  of  '  were 
fo\ind: 


n 

1 

2 

3 

4 

1 

0.280074 

0.528107 

0.752163 

0.956771 

2 

0.280074 

0.528107 

0.752163 

0.957406 

3 

0.280074 

0.528107 

0.752163 

0.957405 
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